Scattering of Phonons by a Vortex in a Superfiuid 
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Recent work gives a transverse force on an isolated moving vortex which is independent of the 
normal fluid velocity, but it is widely believed that the asymmetry of phonon scattering by a vortex 
leads to a transverse force dependent on the relative motion of the normal component and the 
vortex. We show that a widely accepted derivation of the transverse force is in error, and that a 
careful evaluation leads to a much smaller transverse force. We argue that a different approach is 
needed to get the correct expression. 
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I. INTRODUCTION 

Quantized vortices have been an essential part of the 
theory of superfluids for the best part of the last few 
decadeaH. Vortices provide an important mechanism for 
the decay of supercurrents in either superfluids and su- 
perconductors, and their motion is responsible for dissi- 
pative effects. A considerable effort has been devoted to 
the derivation of the equation of motion for vortices. 

It is a widely held concept that for a vortex moving 
with velocity vy through a superfiuid there is a force 
transverse to its velocity in analogy to the Masnus or 
Kutta- Joukowski force in classical hydrodynamicsa. This 
transverse force is given by the usual expression 



p s k q x (vy - v s 



(1) 



where p s is the density of the superfiuid component, Kq is 
a vector parallel to the vortex whose magnitude is equal 
to the quantum of circulation h/m, and v s is the su- 
perfiuid velocity in the vicinity of the vortex. At zero 
temperature the superfiuid density equals the total fluid 
density and this force is all there is to the story. At finite 
temperatures, however, part of the fluid is in the nor- 
mal component and there is the possibility of additional 
transverse forces proportional to the velocity difference 
relative to this normal component, (vy — v„). The ex- 
istence and value of this additional transverse force has 
been quite controversial and very^different expressions for 
it can be found in the literatureoH. n r. 
Recent work by Thouless, Ao and NiuQ and WexlerB, 
unambiguously shows that the transverse force per unit 
length on a single quantized vortex in a superfiuid is given 
by the universal expression in eq. [l], without any addi- 
tional transverse forces. The argument is based on a 
global argument that relies on the equilibration of the ex- 
citations, so it does not deal separately with the change 
in momentum due to phonon-vortex scattering and the 
modification of the phonon-phonon scattering due to the 
vortex motion. In essence the argument is as follows: 



Wexler 's thermodynamic arguments shows that the co- 
efficient of Kq x v s is indeed — p s . This result does not ap- 
pear to be controversial. Thouless, Ao and NiuQ showed 
that the part of the transverse that depends on the vortex 
velocity vy is given by 



(p s K s + p n K n ) X Vy 



(2) 



where p n is the density of the normal fluid component 
and k s , K n are the equilibrium circulations of the super- 
fluid and normal fluid parts far away from the vortex 
core. The circulation of the superfiuid is clearly the quan- 
tum of circulation kq = h/m and the normal, fluid is as- 
sumed irrotational n n = due to its viscosityH. Since the 
coefficients of vy and v s are equal and opposite Galilean 
invariance leaves no room left for any additional trans- 
verse force depending on the normal fluid velocity v„. 

The arguments presented are rather straightforward 
but there is still considerable belief in the existence of 
additional forces. One should then consider the alterna- 
tive calculations more closely and ask what is the origin 
of this controversy. In part, the confusion arises from dif- 
ferent interpretations of the role played by the scattering 
of excitations the vortex. This is discussed next. 



II. SCATTERING OF EXCITATIONS BY A 
VORTEX 

There is no doubt that scattering of excitations by a 
vortex produces a longitudinal force proportional, 
relative velocity of the excitations and the vortexlM 
For either phonons or rotons this force is quadratic in 
the circulation k, and at temperatures low enough that 
phonons become dominant is fifth power in temperature, 
if only phonon-vortex scattering is taken into account and 
the perturbation of the phonon-phonon scattering due to 
the vortex motion is ignored. 

It is widely believed, however, that there is also a 
transverse force due to an asymmetry in this phonon- 
vortex scattering. The details of this are not clear to us 
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from Iordanskii's papersEl, but a derivation by ClearyEl 
which gives a transverse force linear in k and fourth 
power in temperature (like the normal fluid density due 
to phonons) is .-widely quoted, for example by SoninEl. 
Demircan et aZ.til, on the other hand, find that this force 
vanishes. |-..— , 

It has been arguecfllla that the divergence of the scat- 
tering amplitude at small angles invalidates arguments, 
such as those in reference [l4|, based on the Born approx- 
imation for the scattering amplitude. In fact the (9 _1 
divergence of the scattering amplitude for small scatter- 
ing angles 8 gives no trouble with the longitudinal force, 
which depends on the scattering cross section weighted 
with (1 — cos 8), and only a formal divergence for the 
transverse force, where weighting of the differential cross 
section with sm8 leads to a principle part integral and 
the need for a regularization procedure. 

Examination of the Cleary derivation shows that the 
derivation is in error. Whereas Cleary obtained a trans- 
verse force per unit length proportional to 



3 J3 



(v v - v n )n T /% c 



(3) 



a correct evaluation of the formula he used for phonon- 
vortex scattering gives an answer which term by term is 
of order 



{v v ~ v n ) K 3 T 6 /h 5 c 9 



(4) 



Cleary reached eq. ||| by adding and subtracting diver- 
gent series to his expression. Demircan et al. did not 
get eq. ^ because they did not keep enough terms. They 
did not get eq. |^ because it does not follow from Cleary's 
starting-point. 

Let us consider a phonon in a superfiuid with an infi- 
nite rectilinear vortex on the z-axis. We assume a sta- 
tionary vortex for the calculations that follow. Galilean 
invariance can be used later to calculate the force in a 
more general situational and the oscillatory motion of the 
vortex due to the interaction with the phonon is not im- 
portant for the determination of the transverse force (see 
Sec. Ill D in Ref. |J). Phonons interact with the gauge- 
like superfiuid velocity field and with the density varia- 
tions near the vortexpGpre. In the long-wavelength limit 
the former dominates! 2 ] and phonons satisfy the following 
equation of motion 



d 2 i 



W + 2 v s0 (r) • V 



di 
~dt 



(5) 



where ip represents the phonon wavcfunction and 
v so(?") = e0fto/(27rr) is the superfiuid velocity around 
a vortex with circulation kq. Given the symmetry of the 
problem and the simple form of eq. |^, it is possible to cal- 
culate the differential cross section for a phonon-vortex 
event exactly by using a partial-wave decomposition. 



ipi{r)e 



ik z z+U8—iujt 



(6) 



where u = c\k\ 
(for/^0) 



Cy/k 2 + fcf. Each partial-wave satisfies 



fin 

dr 2 



r dr 



(I + kn /2ircf 



in + k 2 r in = 



(7) 



where a— .small term proportional to v 2 has been 
dropped!! The solutions can be obtained in terms of or- 
dinary Bessel functions J Vl [k r f), where v\ = \l + kKo/2i:c\ 
for / 7^ (note that from eq. 0, Vq — but since s-waves 
have considerable overlap with the vortex core, in general 
one should keep this case apart). One should also note 
that this is not a Born approximation solution but rather 
the exact solution of the phonon wave equation. 

Upon consideration of the asymptotic behavior of the 
Bessel function one can readily obtain the scattering am- 
plitude 



f{0) 




J2 e lS ' +U9 sinS l , 



where 



s i = -n) =?ysign(0 , 



(8) 



(9) 



and rj = —kn^/Ac. Note that for a vortex in superfiuid 
helium \rj\ <^ 1 even for the most energetic phonons. 

The series in eq. || are not convergent, since all the 
terms are of the same magnitude, but when we regular- 
ize the series in an obvious way, by writing the phase 
shifts as 

Si = ?/ sign(^) exp(-a|/|) , with < a < 1 , (10) 
we get the scattering amplitude as 



/reg(#) 



2 shir] cos(r] + 8/2) 
irk r sin(#/2) 



(11) 



in the limit a — * 0. Here the phase shift So has been 
taken to be zero. This gives the differential scattering 
cross section as 



■sin 2 r] [cos 2 !] cot 2 (6>/2) 



da _ 2 
dd irk r 

— 2sin?7 cos rj cot(#/2) + sin 2 rf\ 



(12) 



In this form the differential cross section looks innocuous, 
giving a well-defined transport cross section in terms of 
the integral weighted with 1 — cos#, and a transverse 
scattering cross section, in terms of the integral weighted 
with sin 8, that is a principal part integral. As we will see, 
a different limiting process can, however, give a different 
answer for the transverse scattering cross section, since 
the limiting behavior at very small angles can change the 
integral. 

If we calculate the transverse cross section directly 
from eq. H this gives 
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p2tt 

cr± = d9 sin 
Jo 
2 



l/WI 2 (13) 
amOde V e l{l - l ' )e e l( ~ Sl - s ^ sinc5 ; smS v 



1,1' 



— ^sinJ; sin^+i sin(<5; — 8i+i) 



(14) 



where the sums and integral were commuted to obtain 
the last expression. This last step is not trivial given the 
non-uniform convergence of the series. If the expression 
14 is taken literally, then, for 5q = 0, every term 
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in this series vanishes. However, if the large angular mo- 
mentum phase shifts are smoothly taken to zero, as in 
eq. [Io|, then the sum in cq. [l4] can be expressed in terms 
of an integral over <5; : 



8 

kr JS 

2 



Si 



sin 2 Si sin d5. 



(2?/ — sin 2rf) 



i - 
8rf 



sin 2 5i dSi 



(15) 



If the same integral (eq. [H?]) is performed for the expres- 
sion for the differential cross section given in eq. [l2|, then 
we get 



— sin i] sin 2r\ . 



(16) 



These two expressions, [l5] and |l6|, are both of order rf, 
but differ both in sign and magnitude of the coefficient. 
Examination of the regularization procedure shows that 
in the range of angles a « 9 « y/a the even part of the 
scattering amplitude is proportional to rj 2 a/9 2 , while the 
odd part is proportional to rj/6. The squares of each of 
these terms give contributions to the transport cross sec- 
tion which vanishes in the limit a — > 0, but the product 
gives a contribution to the transverse cross section that is 
of order rf . Let us look at this in more detail. We start 
with eqs. |^ and ^ and assume the same regularization 
procedure as above, expanding this time in powers of r\ 
while keeping a finite throughout the calculation: 



/(#) = t/-7- Y^i sin(?7e- Qi ) shxhe 

' 2>1 



V 



-od 



+ 19) (17) 



cosh a — cos 9 



cosh 2a — cos 



From these expressions, one can calculate the differential 
cross section. The square of the first term in f(9) gives 
the even part of the cross section 



sin 



da e 2 r/ 2 
d9 irk r (cosh a — cos ( 



(18) 



which contributes to a finite transport cross section but 
gives a vanishing transverse cross section <r±. The next 
order term is odd in r\ and 9 



da 4 if sin 9 (cos 9 — cosh 2a + sinh a) 
d9 irk r (cosh a — cos 9) (cosh 2a — cos 9) 



(19) 



and does contribute to a^. If one considers the limit 
a — > at this point, the cross section reduces to the 
expression shown in cq. [l^, and the transverse cross sec- 
tion would be given by eq. |l6|. The problem is that a 
series expansion of da /d9 in powers of a yields diver- 
gent contributions to a± coming from small angles. The 
expression above, however, is completely regular for fi- 
nite a. One then calculates the transverse cross section 
first, which yields a regular function of a whose series 
expansion is given by 



3 A 1 



a 
3" 



' 24fc r c 3 



(20) 



in complete agreement with eq. [Lq. The reason for the 
apparent contradiction between eqs. |l5| and |l6| comes 
from important small-angle contributions which were ne- 
glected in the simple-minded approach leading to eq. [If]. 
It is clear, however, that no contribution of order lower 
than if' is present. 

Why does Cleaut's calculation, give a very different an- 
swer? In Cleary'sE£l and Sonin'sB calculations, eq. |lj was 
set equal to 



crj_ 



l]Tsin2(<5 ; 



1+1 



), 



(21) 



which can be obtained from eq. [l4] by adding and 
subtracting the divergent series (l/k r ) ^ sin2<5; and 
(l/fe r ) 2 sin25;_)_x- A result linear in r\ is thus obtained 
from an expression that is proportional to rf . 

The transverse force is usually calculated by inserting 
the transverse cross section (eq. po|) into the phonon cur- 
rent (note that for this calculation we have vy = v s = 0) : 



d 3 k 
J2nf 



a±(k) n(k) (Tick x e 2 



= -(e 2 xv„) 



67r 2 fcsT 



dk <7±(k) 



}A e hcklk B T 
JgHck/ksT _ ]12 



(22) 



Without going into detailsH it is clear that the fc 2 depen- 
dence of <7 _l (eq. p0| ) yields a transverse force proportional 
to T 6 as stated in eq. [|. In contrast, the normal density 
is proportional to T 4 as can be clearly seen from 



• ph 



d3k *L. M 2tt 2 (k B Tf 

Tik n(k) ~ — — „ v„ = p n v n . (23) 



(27T) 



45 Ti 3 c 5 



We have seen that within this framework the result is 
of order rj 3 , but that the coefficient of this term depends 
on the limiting processes which are used. There are other 
reasons that make the result suspect. The first is that 
to the order of eq. [| there are corrections to the sim- 
plest hydrodynamic expression for the phase shifts which 
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might contribute to the same order, but which cannot 
possibly contribute to a lower order in T . The second 
is that that there should be contributions from the s- 
wave phase shift So, which will depend on the details 
of the vortex core structure. It is easy to see from eq. 
[14] that for So 7^ there is a contribution to a± given by 
(4/Av) sin2?7 sin 2 So ~ ko<5q. While usual scattering prob- 
lems in two dimensions have s-wave phase shifts that de- 
pend logarithmically on the momentum k, for this prob- 
lem the boundary conditions are suchE3 that So ~ (k r a) 2 , 
where a is of the order of the vortex core radius. In this 
case the additional contribution to the transverse force 
will be of order 



(v v - v n )n T /h c 



7 A3 



(24) 



Although eq. [| might be canceled by corrections to 
the hydrodynamic theory, there is no chance that eq. |24| 
can cancel exactly, since it is strongly dependent on the 
details of the core. 



III. CONCLUSIONS 

Our analysis has— shown that a correct evaluation 
of Cleary's formula!!! for the transverse scattering of 
phonons by a moving vortex gives a transverse force much 
smaller than the one he obtained, and which seems to be 
widely accepted. However, it is nonzero, and slightly 
sensitive to the details_af the vortex core, which seems 
to contradict our clairrofl that the transverse force on an 
isolated vortex is independent both of the structure of 
the vortex core and of the normal fluid velocity. 

This is a point that requires more detailed study, but 
we believe that the Cleary calculation, involving the sin# 
weighted average of the differential cross section, omits 
the most important effect of the vortex on the phonons. 
This approach assumes that the incoming wave is asymp- 
totically free, but actually for a system of radius R each 
angular momentum mode has an energy shift cSi/nR, 
and so this phase shift leads to a different equilibrium 
population of positive and negative I modes. The shift 
is inversely proportional to R, but the number of modes 
contributing is proportional to R, so this effect is not 
obviously negligibleEa. 

The only way we know to take this into account is 
to use the method which we used in Ref. ||, but with 
more careful examination of the boundary conditions 
than was given in that paper. We are working on 
a careful evaluation of the problem of noninteracting 
phonons in a moving vortex, and our preliminary results 
suggest that the phonons give a negative contribution 
close to —p n Ko x vy, which must be added to the zero- 
temperature result pKo x vy . This term comes from the 
negative net circulation of the phonon momentum around 
the vortex. This results need to be checked carefully, 
since it is in accord with our own preconceptions. 
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